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Abstract 



We have simulated the two-dimensional RP 2 and RP 3 a- models, 
at correlation lengths up to about 220 (resp. 30), using a Wolff- type 
embedding algorithm. We see no evidence of asymptotic scaling. In- 
deed, the data rule out the conventional asymptotic scaling scenario 
at all correlation lengths less than about 10 9 (resp. 10 5 ). Moreover, 
they are consistent with a critical point at (3 ~ 5.70 (resp. 6.96), only 
2% (resp. 5%) beyond the largest (5 at which we ran. Preliminary 
studies of a mixed S N ~ l /RP 1 ^ -1 model (i.e. isovector + isotensor 
action) show a similar behavior when (5t — > oo with (5y fixed < 0.6, 
while they are consistent with conventional asymptotic freedom along 
the lines Pt/Pv fixed < 2. Taken as a whole, the data cast doubt on 
(though they do not completely exclude) the idea that RP 1 ^' 1 and 
S N_1 cr-models lie in the same universality class. 
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The present work grows out of our study of Wolff-type embedding algo- 
rithms for general nonlinear cr-models ]l|, ^ . Our main result was that a 
Wolff-type embedding algorithm can perform well — in the sense of having a 
dynamic critical exponent z -C 2 — only if it is based on an involutive isom- 
etry of the target manifold M whose fixed-point manifold has codimension 
1. Such an involutive isometry exists only if M is a sphere, a real projective 
space, or a discrete quotient of products of such spaces. 

After extensive studies of the models where M is a sphere S"^ -1 (the so- 
called N- vector model), we considered next the models with M = RP N ~ 1 = 
S N ~ 1 /Z2 ||. We were thus led to re-examine the long-standing puzzle con- 
cerning the phase diagram of these models. 

This interest is not purely academic. RP 2 a-models are used as models for 
nematic liquid crystals, where instead of spins one deals with the orientation 
of elongated molecules exhibiting a head-tail (<x — > — cr) symmetry. 

On a Z 2 lattice we shall consider (for the moment) 

- (3H T ({a}) = ^EK-^f W 

(%y) 

(T for "tensor") with <r x G S N ~ l . Clearly is invariant under the Z 2 gauge 
transformations 

a x -> r] x cr x (2) 

with T) x G Z2 = {±1}- It follows that all non-gauge-invariant correlation 
functions vanish identically, for example 

(0Z0*) = if x ± y . (3) 

The fundamental correlation function is thus the isotensor 

G T (x,y) = ((* x -cT y ) 2 )-l/N . (4) 

In the formal continuum limit, the discrete gauge symmetry plays no role. 
Locally the spaces S 1 ^^ 1 and RP 1 ^ -1 are the same, and one thus expects the 
model ([!]) to have the same continuum limit in the isotensor sector as the 
ordinary iV-vector model. 

At iV = 00, the model ([!]) undergoes a first-order transition j|; but this 
transition involves the spontaneous breaking of the Z2 gauge symmetry, and 
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thus cannot persist to finite N. For finite N, some workers have found (by 
high-temperature expansions or Monte Carlo simulations) indications of a 
second-order phase transition, which could be related to a condensation 
of Z 2 -vortices in the field 

An(x) = sign [er x ■ er x+fl ] . (5) 

If there is not a phase transition, the model is presumably asymptotically 
free, and a renormalization-group calculation predicts that for N > 2 the 
correlation length £ and susceptibility x behave as 

C f (/3) ~ C ^~l/(N-2) e 2^/(N-2) (6) 
X af (P) ~ Cx(3 -(N+l)/(N-2) e W/(N-2) (?) 

up to corrections of order 1/(3. 

By using the exact results for the mass gap for the 0(iV)-models and 
the ratios of A-parameters in different schemes (which are extracted from 
first-order perturbation theory), it is even possible to predict what is the 
value of the constant for the exponential correlation length (inverse mass 
gap). In the following we shall make use of the so-called second- moment 
correlation length, which should differ by a few percent from the exponential 
correlation length. 

In Fig. [I] we report the ratio of the measured correlation length to the 
value (0) predicted by asymptotic freedom, as a function of (3 on lattices 
of sizes 32 < L < 512, for both RP 2 and RP 3 . The observed correlation 
length is smaller than the expected value by a factor of about 10 7 (resp. 
10 4 ). The simulations at the larger values of (3 and L clearly show that the 
asymptotic regime has not been reached, as the curves show a clear rise up. 
In any case, our simulations on small lattices prove that if the conventional 
scenario is true, the asymptotic regime cannot be seen on lattices of size less 
than about 10 9 (resp. 10 5 ) — a result which is quite astonishing. 

Let us assume, instead, that there is a phase transition at a finite (3 = (3 C . 
Then we expect by finite-size scaling that 

m = Lfr(\p-p c \ v L) (8) 
X ((3) = L^f x {\f3-(3 C \»L) (9) 

where L is the size of the lattice, and v and 7 are, as usual, the critical 
exponents for the correlation length and the susceptibility. Thus we could 
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try to fix j3 c and v by using the data for £ on different lattice sizes and by 
requiring that all the obtained values of £/L, when plotted as a function of 
(1 — P/PcjL 1 ^, fall on the same curve. Analogously, once (5 C and v have been 
fixed, 7 can be found by imposing @. 

In Fig. H we have plotted the observed values of ^/L as a function of 
(1 — (3/P c )L 1 / u with P c = 5.70 and v = 2.0 (see also ||). The same analysis 
for x fi xes 7 ~ 3.3. Also for RP 3 it is possible to fit the data by supposing 
a /3 C — 6.96 and the same values of v and 7 as for RP 2 . 

The quite large values of the critical exponents may suggest a phase 
transition of the Kosterlitz-Thouless type: in this case the behaviour given 
in (|) and @ should be substituted by 

= L/ e (exp[-A(l -(3/(3 c )- 1/2 \ L) (10) 

X {0) = L^f x (exp[-A(l - (3/(3 c r 1/2 } L) (11) 

where fllCf ) should be used to determine the constants A and f3 c . In Fig. ^ we 
have plotted the observed values of £/L as a function of exp(— A/\l — (3/ (3 c \ l ^ 2 )L 
with (3 C = 5.85 and A = 1.73. The same analysis for \ fixes 7/V ~ 1.7. Also 
for i?P 3 it is possible to fit the data by supposing a (3 C = 7.15 and the values 
A = 1.7 and 7/1/ = 1.65. 
A plot of x/£ 7/l/ 

clS cl function of £/L provides directly an estimate of 
7/z/, which is independent of the assumptions on the scaling behaviour. See 
Fig.|. 

Next we considered the Hamiltonian 
—{3H mixed ({cr }) = 

= Pv X) a * • °» + y X( CT ^ • °"y) 2 ( 12 ) 

(V for "vector") which interpolates between the models with target space 
S 1 ^^ 1 and RP N ~ X . The normalization of the two terms has been chosen 
in order to get (3 = Pv + Pt as the coupling constant in the usual low- 
temperature expansion. 

For iV = 4, we made simulations along the lines with r = Pt/Pv fixed and 
not too large (0.5, 1.0, 2.0). We found that the model behaves qualitatively 
in agreement with the conventional scenario of asymptotic freedom in the 
sense that 
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• the ratio £/£ a ^ is not very far from 1, although it becomes worse and worse 
when r increases and at r = 2.0 it is around 0.4; 

• the estimated values of (3 C which appear when we try the Ansatze (|8|) / (|9]) or 
( |10"D / (|ll|) are relatively far from the actual values of f3 which can be simulated 
effectively. But this gap decreases as r increases; 

• the ratio of the vectorial over the tensorial correlation length seems to ap- 
proach a constant value in the continuum limit; 

• the dynamic behaviour of the Wolff embedding algorithm is in the same 
universality class as for the O(N) model, which means that the dynamic 
critical exponent is almost zero (critical slowing-down is eliminated). 

On the other hand, in the region where fix — > oo with f3y fixed and small 
(< 0.6), the model looks very similar to the RP N ~ X case: the ratio 
is extremely small, and the estimated values of a transition point (3 C are 
very close to the largest values of (3 at which we ran. The vectorial sector 
appears decoupled from the tensorial one: the vectorial correlation length and 
susceptibility appear to approach finite values even when the corresponding 
tensorial quantities start to diverge. Also from the dynamical point of view 
the behaviour is the same as in the RP 1 ^^ 1 case, with a dynamic critical 
exponent for the tensorial susceptibility close to 1 0. 

But we have also seen, in our Monte Carlo runs for N = 3, that there is 
a crossover, near (3y w 0.5, where also the vectorial sector begins to become 
critical. Ising-like fluctuations seem to become important in this regime, so 
that the Wolff embedding algorithm is less efficient. We are currently in- 
vestigating this region: it would be interesting to understand whether the 
apparent transition curve which starts on the axis (5y = stops somewhere, 
and whether the Ising critical point which occurs at (fix — 00 , Pv — 0.44 . . .) 
is isolated or belongs to a phase-separation curve extending to finite fix- The 
phase diagram could appear like in Fig. ||, but are the phase transitions ac- 
tually there? An unconventional scenario has been suggested || in which 
asymptotic freedom does not hold even in the pure vector model. Indeed, 
conjectured correlation inequalities, based on the remark that the couplings 
are both ferromagnetic, would suggest that the line where the tensorial cor- 
relation length diverges — if there is one — cannot stop somewhere in the 
middle of the phase diagram, and must at least intersect every ray with 
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r = Pt/Pv fixed > 0. We are continuing to run our simulation, and hope to 
add other pieces of evidence to clarify these important open questions. 
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Figure 1: Asymptotic freedom test for RP 2 and RP 3 . 



Figure 2: Scaling for the correlation length for RP 2 . 



Figure 3: Scaling of the Kosterlitz-Thouless type for the correlation length 
for RP 2 . 



Figure 4: x/L 1 ' 7 vs. £/L for iV = 3 with < (3 V < 0.5. 



Figure 5: A possible phase diagram. 
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